Parasitic interference in long baseline optical interferometry - Requirements for 

hot Jupiter-like planet detection 

A. Matter 1 , B. Lopez, S. Lagarde 

Laboratoire Fizeau, UMR 6525, UNS - Observatoire de la cote d'azur, BP 4229, F- 06304 Nice 

Cedex 4, France 

W.C. Danchi 

NASA/GSFC, Greenbelt, MD 20771, USA 

S. Robbe-Dubois, R.G. Petrov 

Laboratoire Fizeau, UMR 6525, UNS - Observatoire de la cote d'azur, 06108 Nice cedex 02, France 

and 

R. Navarro 

NOVA-ASTRON, P.O. Box 2, 7990 AA Dwingeloo, The Netherlands 

Abstract 

The observable quantities in optical interferometry, which are the modulus and the 
phase of the complex visibility, may be corrupted by parasitic fringes superimposed 
on the genuine fringe pattern. These fringes are due to an interference phenomenon 
occurring from straylight effects inside an interferometric instrument. We developed 
an analytical approach to better understand this phenomenon when straylight causes 
crosstalk between beams. 

We deduced that the parasitic interference significantly affects the interferometric phase 
and thus the associated observables including the differential phase and the closure 
phase. The amount of parasitic flux coupled to the piston between beams appears to be 
very influential in this degradation. For instance, considering a point-like source and a 
piston ranging from A/500 to A/5 in L band (A = 3.5 /im), a parasitic flux of about 1% 
of the total flux produces a parasitic phase reaching at most one third of the intrinsic 
phase. The piston, which can have different origins (instrumental stability, atmospheric 
perturbations, ...), thus amplifies the effect of parasitic interference. 
According to specifications of piston correction in space or at ground level (respectively 
A/500 ~ 2nm and A/30 ~ lOOnm), the detection of hot Jupiter-like planets, one of 
the most challenging aims for current ground-based interferometers, limits parasitic 



Corresponding author: matter@oca.eu 



- 2 - 



radiation to about 5% of the incident intensity. This was evaluated by considering 
different types of hot Jupiter synthetic spectra. Otherwise, if no fringe tracking is used, 
the detection of a typical hot Jupiter-like system with a solar-like star would admit a 
maximum level of parasitic intensity of 0.01% for piston errors equal to A/15. If the 
fringe tracking specifications are not precisely observed, it thus appears that the allowed 
level of parasitic intensity dramatically decreases and may prevent the detection. In 
parallel, the calibration of the parasitic phase by a reference star, at this accuracy level, 
seems very difficult. Moreover, since parasitic phase is an object-dependent quantity, 
the use of a hypothetical phase abacus, directly giving the parasitic phase from a given 
parasitic flux level, is also impossible. Some instrumental solutions, implemented at the 
instrument design stage for limiting or preventing this parasitic interference, appears to 
be crucial and are presented in this paper. 

Subject headings: parasitic interference, mid-infrared interferometry, phase, hot Jupiter 



1. Introduction 

Stellar optical interferometry has substantially evolved in terms of instr ument and op eration 



since Fizeau's idea to use this technique for measuring the diameter of stars (|Fizeaulll868l ). Since 



Michelson's inte rferometer, where separated mirrors were fixed on the same telescope mounting 



(lMichelsorjll920h . in terferometry with independent telescopes has allowed access to very long base- 
lines ( Labeyriel 1975 ) and consequently high angular resolution. In this case, which corresponds to 



most current interferometers, the beams coming from each telescope, are carried through tunnels 
up to a combining device. The beams are often reduced in size for practical reasons, and the cur- 
rent instruments are characterized by optical modules performing various functions such as spatial 
filtering, spectral band separation, and spectral resolution. The transport of these beams through 
multiple optical modules can be problematic. 

A parasitic interference may occur because of diffraction effects associated with beam propagation 
along finite size optics. This diffraction can then produce a cross-talk between beams. The result 
is the superposition on the detector of several fringes systems having different phases. Therefore 
the interferometric signal of the independent beams is affected. While the intrinsic interferometric 
fringe pattern corresponds to Fizeau-like fringes, the parasitic interference induces two further com- 
ponents to the signal: Fizeau-like mirror fringes and Young-like fringes. The latter contribution is 
independent of the object position in the sky. We will develop this description in the article. 
This interference corrupts the value of interferometric observables : modulus and phase of the 
complex visibility. The understanding of this phenomenon is of general interest for optical long 
baseline interferometers which combine multiple beams and often require some "compactness" in 
the opto- mechanical elements. 

In Section 2, a simple formalism is developed, explicitly showing the different contributions of the 
parasitic interference. In Section 3, the resulting interferometric observables are written for the 
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general case of an extended source. Then the theoretical cases of an unresolved source and of a 
stellar system with a hot Jupiter are highlighted. Sections 4 and 5 quantify the impact of the 
parasitic interference on the theoretical phase signal of these types of sources. The parameters 
involved in this quantitative study are the factor of parasitic flux between beams, the piston, and 
the photometric imbalance between interferometric arms. The impact of the parasitic interference 
on hot Jupiter detection is evaluated by using different synthetic spectra of such extrasolar planets. 
Requirements on the straylight level limits are also given. 



2. Formalism of the problem 
2.1. Overview 

In general, long-baseline stellar interferometry consists of sampling an incident wave packet by 
means of telescopes at different locations. These coherent waves are combined and the resulting 
interference pattern is extracted in order to measure the complex degree of coherence of the radiation 
field. Then the brightness distribution of the source can be recovered thanks to the Van Cittert and 
Zernike theorem. This complex degree of coherence corresponds to the covariance of the electric 
fields collected by the telescopes. If a parasitic interference occurs between the collection and 
the recombination steps, the 'intrinsic' coherence between beams and consequently the resulting 
interferometric observables will be perturbed. 

To our knowledge, the issue of parasitic interference has never been formalized except for the 
di gital wave-front meas uring interferometry technique. This well-established technique, described 



m 



Bruning et al.l ()1974l ). allows testing of optical surfaces and lenses, and measurement of wave- 



front deviations in the A/100 range. Some systematic error sources such as 'extraneo us' fringes, that 



i s par asitic interference, were examined in the same paper and more extensively in lSchwider et al 



(j 19831 ). In our case, the theoretical description of this perturbation is hereafter detailed in the 



framework of a simple two-telescope interferometer. 



2.2. Interferometric framework 

To create a model of parasitic interference, we use a two-telescope interferometer characterized 
by a multi-axial scheme and an image plane recombination. Fig. 1 gives an illustration of such 
an interferometer observing an unresolved astrophysical so urce. For a gener al description of the 



output response of a single baseline interferometer, see also lElias et al.l (|2007l ). 
First we write the complex amplitudes collected by both telescopes, i.e. the two samples of the 
wavefront, that we respectively note ipi and ip2- We multiply each of them by a real transmission 
factor noted ti and t2, which represents the transmission of electrical fields through the instrument. 
Let us define e to be the main fraction of ip\ propagating towards the path 1, and e' to be the small 
fraction propagating towards the path 2. We assume the same parasitic effect for rp2. Fig. 2 
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shows two examples of cross-talk between beams. A possible imbalance of parasitic flux between 
paths, not presented in this paper, was considered in the study but appeared to be a second order 
parameter. All the cross-talk occuring inside the instrument produces a resulting parasitized pattern 
in a conjugate of the pupil plane, located just before the recombination on the detector. This is 
described in the following equation giving the complex amplitude preceding the recombination : 

^pu P (x, y) = [eti^i + e't 2 ip2\P{x -^,y) + [et 2 ip2 + ^hipi]P(x + -, y). (1) 

Here x and y are the coordinates in the pupil plane, b is the distance between the pupils of 
both telescopes, which are reduced in size at the entrance of the interferometer, D being the 
diameter of these pupils. Since in a Michelson configuration the collection and correlation steps 
are separated, b, the 'reduced' baseline which forms the fringe pattern, is different from B, the real 

interferometric baseline which samples the incident wavefront. P(x, y) = II( V D — ), where H(x, y) 
is the transmission function of a circular aperture with a uniform transmission of 1 inside and 
outside. From Eq.(l), we discuss how parasitic fringes are formed on the detector. 
As we will see later, the complex degree of coherence between both parasitized beams will become 
< [eiiV>i + e^V^] [^2^2 + e'tiipi]* >, instead of the 'intrinsic' coherence term < V'iV'2 >• 



2.3. Interference and formation of parasitic fringes 

In order to describe the formation of parasitic fringes, we consider the most general case of 
the observation of an extended source. The vector en is the angular coordinate in the plane of the 
sky. The telescopes are located at 71 and 72 , 7 being the coordinate in the plane containing the 
telescopes and counted in units of wavelength (7 = j) (see Fig. 1). ^ = 71 — 72 with B the 
interferometer baseline. The electric field emitted by each point of the extended source, located at 
a, is represented by its amplitude A(a) and its phase 4>(a,t) : 

Y> s (a) = A(a)e^ + ^ a '*)W. (2a) 

In the plane of the telescopes, the phase shift of the wavefront emitted by each point of the source 
and measured on the 7 position is A<E> = — 2na ■ 7. Therefore, when considering the contributions 
of all the emitting points of the source, each telescope observes a packet of wavefronts (respectively 
Y>1 and 1P2) ■ 

Y>i =jA(a)e i ' u '^( a ' i )'e ,2mil (iQ 1 (2b) 
Y> 2 = /A(a)e i ^ + ^ a '*)] e i27ra ^ 2 da. (2c) 

W can now write the expression of the complex amplitude in the detector plane by performing the 
Fourier transform of Eq.(l) with respect to x and y. In the detector plane, where f3 and r] are the 
conjugate angular coordinates linked to x and y, and P(f3,rj) is the pupil diffraction function, we 
obtain : 

W((3, rf) = (etiVi + e'*2^ 2 )P(/?, V)e^ + (et 2 ^2 + e'iiViTO, v)e^ = X, + X 2 . (3) 
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+ and (et2'02 + e'*iV ; i) are n °t affected by the Fourier transform since ipi and ^2 can be 

considered, in a good approximation, as constant over the aperture telescope area. This assumption 
is equivalent to the requirement that the source is not resolved by the telescopes themselves or that 
the field of view is small. In Eq.(3), X\ = (etitpi + e't2ip2)P{P,v) e x and X2 = (^2^2 + 
e ti^i)P(/3,r])e A represent the Fourier transform of the parasitized complex amplitudes of each 
beam. On the detector, we observe the following intensity pattern : 

I(P,ri) =< |Xi| 2 > + < \X 2 \ 2 > +2Re < X X X* 2 > . (4a) 

The photometric terms, < |Xi| 2 > and < \X2\ 2 >, are equal to: 

< |Xi| 2 > = P 2 (/3,r/) < |etiV>i + e't 2 ip 2 \ 2 > 

= P 2 (/3,r/)[e 2 t 2 < IV'iI 2 > +e'H 2 2 < |V> 2 | 2 > +2ee't 1 t 2 Re < fai/,* >} 

= P 2 (/3,r/)[(e 2 i 2 + e /2 t 2 )6(0)] + 2ee't 1 t 2 Re(6( 11 - 72)), (4b) 

and 

<|X 2 | 2 > = P 2 (/3, V )[(eH 2 + e' 2 t 2 )6(0)] + 2ee't 1 t 2 Re(6('y 2 -j 1 ))]. (4c) 

Here 0(a) = A 2 (a) is the intensity distribution of the source. The Fourier transform of the 
intensity distribution of the source at the spatial frequency (71 — 72), also called the complex degree 
of mutual coherence, is : 6(71 — 72) =< V'iV'2 >= / 0(a)e~ t2na '^' yi ~^da. In this covariance 
operation, one of the integrals disappears because astrophysical sources are spatially incoherent 
when time averaged. We also note that O(O) =< \ip\\ 2 >=< |V>2| 2 >= f 0(a)da is the total 
energy radiated by the source and collected by each telescope. 

The term 0(a) is a real function, so its Fourier transform, noted 6(71 — 72) = pi2e** 12 , is a 
complex function with an even real part and an odd imaginary part. Finally, the addition of both 
photometric terms gives : 

< |Xi| 2 > + < \X 2 \ 2 >= P 2 (M[(4 + *l)(e 2 + e /2 )6(0) + 4ee'M2Pi2 cos($i 2 )]. (4d) 

In order to highlight the different fringe patterns appearing in the interferogram, let us explicitly 
show the real part of the correlation term, < XiXg >, containing the coherent flux and the spatial 
modulation : 

Re < X1X2 > = Rc(< P 2 (/3,?7)[(etiVi + e'i 2 V> 2 )e^^][(e* 2 V> 2 + etnpl)e^^] >) 

= P 2 (/3,77)Rc[[e e '(t 2 +tl)6(0) + e 2 t 1 t 2 6( 7 i -7 2 ) + e' 2 tit 2 6*(7i - 72 )]e- 2 '' r ^] (4e) 

= P 2 (/3, r?)[e 2 tii 2 pi2 cos(^M + $ 12 ) + ee '(fi + ^)6(0) cos(^) + e' 2 iifepi 2 cos(^ - <E- 12 )]. 

A A A 

In this expression, each cosine modulation factor has a physical meaning. The first modulation in 
cos(^6/3 + $ 12 ) is the intrinsic fringe pattern. The second modulation in cos( 2 ^6/5) corresponds 
to Young fringes created by the interference between each of the beams and their corresponding 
diffracted part. The position of these fringes is fixed and does not depend on the object position in 
the sky. The last modulation in cos (^6/3 — $12) is due to the interference between the diffracted 
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part of both beams. Mirror-like fringes are thus created with an opposite phase with regard to the 
intrinsic fringe pattern. Fig. 1 shows these different fringe patterns in the case of an unresolved 
source. Each pattern has a different amplitude depending on the nature of the object and the 
amount of parasitic flux. 



3. Highlight of the interferometric observables 

3.1. Resolved source 

After describing the different fringe patterns due to the parasitic interference, we highlight the 
resulting interferometric observables that have been degraded. From the general correlation term 
< XiXg >, we extract the resulting parasitized coherent flux p cxt and parasitized phase Xext : 

< XxX^ > = P 2 (/3, r))[ee'(tl + 4)6(0) + e 2 t 1 t 2 6( 7l - 72) + e' 2 M 2 6( 72 - 7 i)]e" 2 ^ 

= P 2 (/3,77)[ Pcxt e *«*] e" 2 - b f = P 2 (/3, r,) Pext e -^<-^\ (5a) 

The resulting interferogram in the detector plane is : 

J(/3, 7?) = P 2 (/3, r,)[(tl + t 2 2 )(e 2 + e' 2 )6(0) + Aee'hhp^ cos($ 12 )] + 2P 2 (/3, r,)p cxt cos(2^ - x «t) 

= Jo,ext[l + V cxt cos(2tt6^ - Xext)], (5b) 

where 

/o.cxt =P 2 ( / 3,77)[(t?+tl)(e 2 + e' 2 )6(0)+4ee^ 1 ^p 12 cos($i 2 )], (5c) 



2P 2 (/3,77)p oxt 2 J [e"(l + t 2 2 )6(0) + t 12 (l + e"2)p 12 cos($i 2 )] 2 + [t 12 (l - e" 2 )pi 2 sin($ 12 )]2 

Vext = ; = — , (5d) 

^o.ext (1 + e" 2 )(l + t 2 12 )O(0) + 4e"t 12 pi2 cos($i 2 ) 

. r t 12 (l-e" 2 )sin($i 2 ) 

Xoxt = arctan — . (5e) 

e"(l + *? 2 )O(0) + ii 2 (l + ^" 2 )pi2 cos($i 2 ) 

ti2 = |j is the transmission ratio between both arms, e" = ^ represents the percentage of parasitic 
contribution, evaluated with respect to e which is the fraction of the electric field passing through 
the right path, e" 2 = (^) 2 is the equivalent ratio in terms of flux (or intensity), later called 'par- 
asitic flux factor'. io,ext is the photometric factor containing the parasitic contribution. Similarly, 
V cxt and Xext are the parasitized visibility and phase. As a result, the parasitic phase, representing 
the parasitic contribution added to the intrinsic phase $12, is Xext — ^12 • 

To summarize, the photometry, visibility and phase of a resolved source are degraded. This degra- 
dation depends on the complex degree of coherence of the source expressed by p\i and $12. The 
coherent flux p cxt can not be calibrated by using the parasitized photometry even without any 
photometric imbalance (ti 2 = 1). It also appears that the phase can not be calibrated. 
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3.2. Unresolved source 

In this subsection let us consider the parasitic interference effect while observing a point-like 
source. The unresolved object is described by its angular position in the sky noted o:p i n t, with 
respect to the line of sight, and its monochromatic flux noted I po int(A). The chromatic brightness 
distribution can be represented in the same way as in the resolved source case, such that only one 
tilted wavefront originates from the source : 

Opoint(Q:, A) = /point (A)5(a - Qpoint)- (6a) 

Opoint(^ 5 A) is the Fourier transform of the brightness distribution, in the plane of spatial frequency 
i*j(u, v) covered by a single-baseline interferometer. The u-axis is chosen along the interferometric 
baseline defined above and takes the value u = -v-. Therefore we have : 

OpointK A) = I P oi„t(A)e i27ru - a p^. (6b) 

The fringe phase is the argument of OpointC 1 ", A) : 

^point(A) = 27TU • ttpoint- (6c) 

Normally, $ p0 mt(A) should be defined with an additionnal unknown constant. The reason is that 
the measurement of the interferometric phase is affected by an ambiguity on the zero-delay point, 
that is the origin of the fringes, and by an unknown integer number of 2ir phase rotation. Therefore 
the phase can only be addressed in a relative point of view, either between the source and a reference 
object, or between several spectral channels. The latter constitutes the differential approach which 
allows to remove the phase ambiguity thanks to a colour difference with an appropriate reference 
wavelength. < I ) point(A) is considered to be such an unambiguous differential phase, where the phase 
reference has been taken as the phase averaged on all the spectral channels of the observing band. 
Next the coherent flux is the amplitude of O po i nt (u, A) : 

Ppoint(A) = Ipoint(A). (6d) 

The parasitized photometry, visibility, and phase of an unresolved source are deduced by replacing 
$12, P12 and O(0) in Eqs.QQQ Replacing $i 2 , p\i and 0(0) by $ p0 mt(A), p po int(A) and I p0 mt(A) in 
Eq.(5c), Eq.(5d) and EcL(5e), allows us to deduce the parasitized photometry, visibility and phase 
of an unresolved sourcqj : 

Jo, P omt(A) = JpomtP 2 (/3, ri)[(e 2 + e 2 ){t\ + t%) + Uttxti cos($ P omt)], (6e) 



T/ m , - VQ- + t" 2 ) 2 tj 2 - igtfg Sin^(c& p0illt ) + C »2(l + tj 2 Y + 2tl 2 (l + 6" 2 )6"(1 + t\ 2 ) COs($p 0mt ; 

V »**W- 2 (t? 2 + l)( e ''2 + l) + 4e''t 12 cos(4. poi ^ ' m 



, / sin($ P oint)(l - e" 2 )ti2 s lR 

Ypoint(A) = arctanl ^— - — ; ' — ; ?r-r). de 



The A dependency of $p int(A), p po int(A) and I po int(A) has been removed in Eq.(6e), Eq.(6f), and Eq.(6g), for 
lightening the notations. 
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Here 7i — 72 = x an d ot are considered to be colinear so we have : $ po int(A) = 2n ^ a . Regardless 
of the angular position of the star in the sky, other sources can contribute to the phase term, 
< &point(A), such as instrumental and/or atmospheric perturbations. As shown in Fig. 1, the total 
piston noted L and resulting from these different causes, is first corrected by the delay lines with an 
OPD model depending only on the angular distance between the source and the zenith. Then the 
remaining piston, noted by A in Fig. 1, may be corrected by a fringe tracker which commands the 
delay lines to adjust to a value of A with an accuracy defined by the fringe tracking specifications. 
This accuracy is noted d and represents the residual piston that will affect the fringes within the 
signal envelope. Therefore the generalized expression of the phase term will be : &pointW = 
Similarly to the resolved case, we consider Xpoint(^) — < I ) point(A) as the parasitic phase. 
According to Eq.(6f), Vpomt = 1 if ti2=l, and depends neither on e" and e" 2 , nor on the piston. In 
the case of an unresolved source, the degradation due to the parasitic light can thus be corrected 
if we divide the parasitized coherent flux by the parasitized photometry. Though, this correction 
is effective only if the photometric imbalance ti2 is close or equal to 1. 

3.3. Hot Jupiter-like source 

The direct detection of Earth- like extrasolar planets and the determination of their atmospheric 
features by means of stellar interferometry is a very challenging objective. Nowadays the observation 
of hot Jupiter-like planets constitutes an intermediate step in terms of difficulty. These planets are 
giant gaseous planets similar in size to Jupiter except that they orbit very close to their star. This 
proximity provides planet temperatures warm enough to allow significant infrared excess. The 
planetary flux can be decomposed in this way : 

/pi (A) = /pl,rcfl(A) + /pl,abs(A) + /pl,int(A), (7) 

where / p i, re fl(A) is the component of the stellar light reflected by the planet, f p i )a bs(^) the absorbed 
and thermally re-emitted component, and f p i,int(A) the intrinsic contribution of the planet. For 
non- irradiated objects, the main contribution would be limited to the intrinsic flux f p i : - m t = crT^ t , 
where / p i,int is now the bolometric intrinsic flux and T; n t is the intrinsic temperature of the planet 
(Tjnt = T e ff in this case). For an irradiated planet, T e g is not relevant because it is difficult to 
separate photons which are thermally radiated by the planet from those of the star itself reflected 
by the planet. Another temperature, describing the equlibrium state of the planet's dayside and 
separated from the reflected component, is thus defined : dT c 4 q = c7] 4 t + / p i )a bs) where / p i ja bs is a 
bolometric flux. 

In the near and mid-infrared domain, the major contribution of the flux of hot Jupiter-like exoplan- 
ets is the absorbed, then re-emitted, component at the thermal equilibrium. In purely quantitative 
terms, the equilibrium temperature of such a planet approximately ranges from 500 to 1500 K. For 
the closest objects, the corresponding flux ratio between the planet and the star is about 10~ 5 to 
10~ 4 in J band, 10~ 4 to 10~ 3 in K band, and several 10 -3 in N band. In addition these planets 
exhibit several characteristics that we will detail in Section 4.1. 
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In this context let us consider a stellar system with a planet. We respectively note I* (A) and I p i(A) 
the monochromatic flux of the two components separated by an angular distance a p \. The angular 
position of the center of the star is a*. We first assume that the spatial location and extension of 
each component is independent of the wavelength. The star appears as a disc of constant angu- 
lar diameter D*, defined by the normalized function n(jj-) equal to one inside the disc and zero 
outside. The planet appears as a point-like source defined by S(a). Both spatial distributions are 
respectively weighted by I* (A) and I p i(A). The chromatic brightness distribution of this system is 

O syst (a, A) = I, (A) n(^-^) + / pl (A) 6(a - (a* + a pl )). (8a) 

The Fourier transform of the brightness distribution of the system in the plane of spatial frequency 
u;(u,v), covered by a single-baseline interferometer, is : 

6 syst (u, A) = ^(u)/*(AK 2 ™- Q * + I pl (\)e i2 ™< a *+ a ^ . (8b) 

The u-axis is chosen along the interferometric baseline defined above, and takes the value u = 
y =7i — 72- The intrinsic visibility of the star is k*(u) = C*(u) e J(I>star = 2 ^p *"^ , with J\ 
the Bessel function of the first kind of order one and D* the stellar angular diameter. Since we 
are in a regime where the star is partly resolved or even unresolved by the interferometer, the 
stellar complex visibility is always real and positive, with k*(u) = C*(u) (or identically C*(A)) and 
*£star = 0. Otherwise, if the star was resolved, $ s t ar could be equal to ir. 

In principle, visibility and phase carry complementary informations, and both might be useful for 
detecting and characterizing extrasolar planets. Nevertheless, considering the curre nt possibilities 



of ca libration of both quantities in terms of instrumental and atmospheric effects (jVannier et al 



2006), we prefer focussing our attention on the differential phase observable. The fringe phase is 



the argument of d syst (u, A) E|: 

, m CM sin($«) + J p i sin($, + <frpi) 

•PsysttA) = arctan— — — — (8c) 

y v ; C*/* cos(<I>*) + I p icos($* + $ p i) ! v ' 

with ^p^A) = 2iru ■ ctpi and ( I > *(A) = 2ttu ■ ct*. Identically to the unresolved case, we chose to 
define $ sys t(A) as a zero- mean differential phase in order to remove the ambiguity related to the 
measurement of the phase fringes. We also generalize the expression of ^(A) by taking into account 
other sources like instrumental instabilities or atmospheric perturbations that may contribute to 
this phase term. We thus have : $*(A) = x^*' wriere ^* i s the same residual piston shown in Fig. 
1 and described in the unresolved case. 
The coherent flux is the modulus of O sys t(u, A): 



Psyst (A) = Jem + A + 2CJJ pl cos($ p i). (8d) 



2 The A dependency of $ P i(A), $*(A) / P i(A), I, (A) and C*(A) has been removed in Eq.(8c), Eq.(8d), and Eq.(8e), 
for lightening the notations. 
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Finally, replacing $12 and p\2 by <£ syst (A) and p S yst(A) in Eq.(5e), allows us to deduce the parasitized 
phase of a stellar system with a hot Jupiter, Xsyst(A) : 

_ f 12 (l-0[a/. S in(^)+Jp lS in(^+$ pl )] 
XsystlA; arctan ^ + ^ )(/ ^ + ^ + + e „ 2)[c ^ cog($+) + ^ + ^ • l»ej 

The parasitized phase, Xsyst(A), contains the intrinsic and parasitic contributions of the phase. 
We thus consider Xsyst(A) — ^systC-X) as the parasitic phase since it represents the parasitic phase 
contribution added to the intrinsic phase of the source, <& syst (A). This parasitic term clearly depends 
on <&*(A) and <J> p i(A), which appears both in <£ syst (A) and Xsyst(A) and does not disappear during 
the subtraction. The parasitic phase contribution appears to be object-dependent. 



4. Quantitative study of the 'unresolved source' case 

In this Section, we evaluate the impact of parasitic fringes for a tilted wavefront originating from 
an unresolved source having a projected angular distance to the zenith different from zero. This 
very common situation can imply an uncorrected residual path difference between both telescopes. 
The principle of astrometric measurement of the position of an off-centre star via the phase of 
fringes may be affected by this situation. 

The L band is chosen in this study since it corresponds to the most favourable ba nd for detection of 



several astrophysical sources including hot Jupiter- like exoplanets (jVannierl 120031 ). The parameters 
involved here are the parasitic flux factor e" 2 , which describes the percentage of parasitic flux, and 
the piston 5. This piston will be hereafter expressed in fraction of the L band central wavelength, 
namely 3.5 //m. 

A typical baseline of 100 meters, and a range of piston values, from A/500 to A/5, are considered for 
calculating the phase signal of the source $ po int(A). In this range, two piston values are especially 
of interest : A/500 and A/30; they approximately correspond to the typical specifications of piston 
correction achieved by a fringe tracking device in the case of space and groud-based interferometers 
(respectively 2 and 100 nm). In parallel, for each piston value, a range of parasitic flux factors, 
from to 10 %, is used for evaluating the parasitized phase signal Xpoint(A). 

Finally we calculate a $ = max($p 0in t(A))-min($p 0int (A)), anda par * = max(x P omt(A)-$point(A))- 
min(x P oint(A) — $ P oint(A)). These quantities represent the amplitude, on the extent of L band, of the 
intrinsic and parasitic phases. They are used for comparing the overall amplitude of the intrinsic 
and parasitic phases in L band, instead of comparing both phases at all the wavelengths. As the 
value of the parasitic phase depends on two independent parameters and varies with respect to the 
wavelength, this overall amplitude gives a typical value of the expected measurable phase signal in 
L band. Then this typical value can be plotted in a readable 2D representation as a function of 
the parasitic flux level. In Fig. 3, we plotted these two amplitudes in function of e" 2 , considering 
different values of the piston parameter. We can notice that for all the piston values, the amplitude 
of the parasitic phase, a par *, is a smooth crescent function of the parasitic flux factor. Regarding 
the evolution of the intrinsic and parasitic phase amplitudes as a function of the piston value, it 
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appears that q par* always lies below 10 ^ radians while the amplitude of the phase signal of the 
star, a$, increases up to 0.5 radians. 

Several conclusions can be drawn from this, depending on the context and the objectives of the 
observation : 

• First we place in a classical interferometric context where the piston error has to be corrected 
before the addition of frames. A reasonable parasitic flux contribution of about 1%, and an 
amount of uncorrected piston within the specifications of fringe tracking (A/500 and A/30 in 
L band) are assumed. In this case, the additional parasitic phase amplitude always lies below 
the intrinsic signal, and ranges from one third to one fifth of this intrinsic phase amplitude. 

• Second we consider an 'astrometric' approach, where the aim is to measure the phase signal 
created by an off-centre star, with respect to an on-axis one. It appears here that the mea- 
surement of the intrinsic phase amplitude is not degraded by the additional parasitic phase 
contribution especially for large piston values. The ratio between the intrinsic and parasitic 
phase amplitudes ranges from 3 to 13. 



5. Hot Jupiter signal and magnitude of parasitic fringes in L band 

5.1. Characterization of planetary fluxes 

The fine study of the atmosphere of hot Jupiter-like exoplanets requires a characterization 
of their spectra since they may show significant differences from black body models. Contrary 
to the giant planets of our solar system, hot Jupiters are subject to an extremely intense stellar 
irradiation. Consequently, the temperature-pressure profile of the atmosphere is modified and a 
radiative area is go ing to develop and govern the cooling and contraction of the planet interior 



(jGuillot et al.lll996l ). These structural changes of the atmosphere, compared to an isolated planet, 
strongly depend on how the stellar irradiation penetrates the planet interior. The albedo and the 
sources of opacity, related to the presence of different types of dust, are two important parameters 
allowing us to charact erize and classify the theoretical spectra of hot Jupiter-like exoplanets. 



Sudarsky et al.1 (|2000i ) distinguished five categories of theoretical spectra classified according to a 
range of effective temperatures, from T<150K to T<1500K. The albedo of the objects belonging 
to these classes, exhibits similar features directly related to the chemical composition of the planet. 
In our case, the main criterion we use is the presence or absence of dust in the high layers of 
the atmosphere. For example, in a 'condensed' atmosphere model, dust has been settled down 
by the gravitational field or by the precipitations linked to the condensation and the formation 
of clouds. From a global point of view, this hypothesis is presumably more realistic for substella r 



objects with low temperatures, including the irradiated giant planets (jChabrier and Baraf fc 2000). 
Therefore, for modelling and estimating the variability of the parasitic phase produced during the 
observation of a hot Jupiter- like extrasolar planet, we use eight "condensed" synthetic spectra 
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describing the planetary flux f p i(X). They have been extracted from iBarman et al.l (|200ll ). Each 
spectrum corresponds to a given angular separation and intrinsic temperature Ti nt of the planet, 
which is irradiated either by a G2 or a M5 star respectively located 15 and 5 AU from the Sun. 
The flux of these stars is a simple black body law with temperatures of 5600 and 3000 K. Table 1 
summarizes the values of the different parameters discriminating each synthetic spectrum. 



5.2. Parasitic phase estimation 

In this Section, we estimate the impact of parasitic fringes in the context of the observation of 
different 'synthetic' stellar systems with a hot-Jupiter. According to Sub-Section 3.3, the intrinsic 
phase $ sys t(A), and the parasitized phase Xsyst(A), depend on the flux ratio between the planet and 
the star, the separation p between both components, and the interferometric baseline B; B being 
considered to be colinear with p, and equal to 100m. $ syst (A) and Xsyst(A) are thus calculated for 
each spectrum. The ranges of values taken for the piston <5* and the parasitic flux factor e" 2 are 
the same as those of the 'unresolved source' case. 

As detailed in Fig. 1, <5* can be seen as the remaining piston uncorrected by the fringe tracking 
device. This residual piston creates a linear-like phase term, added to the intrinsic phase signature 
from the planet. Consequently, a further parasitic phase term will be added to the parasitic phase 
term due to the planet phase signal. It is therefore important to place this issue in the framework of 
a classical interferometric measurement since different contributions from the atmosphere and the 
instrument can affect the measurement of the interferometric phase and especially the first-order 
linear term via the residual piston (or piston error) 5*. Basically without any phase reference, the 
first-order linear term of the astrophysical signal is lost during the data reduction step estimating 
and removing this residual piston. This is the case with most current interferometers. Finally, the 
total phase, including the object phase and the parasitic contribution, contains only the higher 
order terms. Henceforth the effect of the parasitic interference has to be examined after removing 
the linear component of <£ syst (A) and Xsyst(A), which consequently become 'unpistoned' phases. 
For a clear illustration of this point, we plotted in Fig. 4 both 'unpistoned' quantities, <3? sys t(A) and 
Xsyst(A)- < I ) syst(A), with respect to the wavelength, using the very smooth number 7 spectrum (see 
Table 1). Xsyst(A) has been calculated for two values of piston (A/500 and A/30) corresponding to 
the fringe tracking specifications, and a quite large parasitic flux factor of 10%. First we can note a 
great similarity between both differential phases, highlighting the fact that the additional parasitic 
phase is an object-dependent quantity. It is quite low and anyway below the intrinsic term in both 
cases. However when the beams undergo a greater piston (here A/15), the parasitic signal exceeds 
the intrinsic one and the phase signature from the planet is lost. 

Identically to the 'unresolved source' case, we consider the amplitude in L band of the phase sig- 
nature from the planet, a syst = max($ syst (A)) — min(<I> syst (A)), and of the corresponding parasitic 
phase, a p ar = max(x sy st(A) — $ sys t(A)) — min(x syst (A) — $ syst (A)). These quantities are calculated 
for a large range of <5* values and are represented with respect to the parasitic flux factor e" 2 in 
Fig. 5. The A/500 and A/30 cases give a very similar parasitic phase amplitude that can only be 
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distinguished for the spectra 5, 6 and 7. Here a detection is considered to be possible only if, for 
given values of (5* and e" 2 , the parasitic phase amplitude lies below the horizontal solid line repre- 
senting the amplitude of the planet signal. Following this criteria, it appears that if typical fringe 
specifications are observed, the planet signal always lies above the parasitic phase amplitude; a hot 
Jupiter detection, with an intrinsic phase amplitude exceeding by a factor of three the parasitic 
one, is thus achieved with a parasitic intensity equal to 5% of the total intensity. 
However, when no fringe tracking is used, the expected constraint on the parasitic flux level appears 
to be much stronger. For instance, let us consider the case of a solar-like star (G2 type) and a 500K 
planet (spectra 5, 6). If the piston error is greater or equal to A/15, a detection, with a planet 
phase amplitude standing at a factor three above the parasitic one, is achieved with a tolerance of 
not more than 0.01% on the parasitic flux. On the contrary, the spectra 3 and 4 constitute very 
favourable cases. If the piston error lies below A/10, a similar detection level admits a tolerance 
of 5% on the parasitic flux. Therefore it appears that as soon as the fringe tracking specifications 
are relaxed or not precisely observed, the parasitic phase signal dramatically increases. In this 
case, planetary detection may be prevented if the parasitic flux level is not accurately monitored. 
Nevertheless in some cases where the flux ratio between the planet and the star is sufficiently large 
(« 10 -2 to 10 -2 ), this assertion has to be moderated given that a non-respect of fringe tracking 
specifications is less influential. In those cases, the tolerance on the parasitic flux level is similar to 
the one achievable with a fringe tracker. 



6. Discussion 

In this work we have described the phenomenon of parasitic interference occuring inside an 
interferometric device. This interference degrades the modulus and the phase of the complex 
visibility. Two further components are added to the intrinsic fringe pattern, the former assimilated 
to Young-like fringes and the latter assimilated to mirror-like fringes with respect to the intrinsic 
pattern. The phase and the amplitude of these 'parasitic' fringes depend on the piston between 
beams, including the piston due to the object position. In an 'object-image' approach, it would 
mean that the point spread function of the interferometer is no longer invariant by translation and 
that the object-image relation is anyway destroyed. 

In quantitative terms, we have shown that the parasitic phase is very sensitive to a piston between 
the interferometric beams. In fact a perfect unresolved source, observed at a projected angular 
distance from the zenith equal to zero, would not produce any parasitic phase. On the contrary, for 
an unresolved source undergoing an additional piston due to the atmosphere and the instrument, a 
parasitic flux factor of 1% would create a parasitic phase reaching at most one third of the intrinsic 
phase amplitude. 

This effect is not negligible in the perspective of hot Jupiter-like planet observations. In this paper 
the feasibility study has been done in L band. We considered synthetic hot Jupiter spectra providing 
different flux ratios between the planet and the star, and different values of piston error within the 
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specifications expected for fringe tracking devices. A hot Jupiter detection, with a planet phase 
signal three times larger than the parasitic phase amplitude, is possible if the parasitic flux reaches 
at most 5% of the total incident flux. Without any fringe tracking device, only hot-Jupiter like 
planetary systems with a quite large flux ratio between the planet and the star, for instance a dM5 
star with a hot planet at 1000K, would admit such a tolerance of 5%; however the values of piston 
error have to be equal or smaller than A/10 in this case. For less advantageous planetary systems, 
the maximum admitted tolerance decreases up to 0.01% for residual pistons equal to A/15. 
To summarize, we can see that the detection of astrophysical objects providing weak signatures in 
the interferometric phase, such as hot Jupiter-like extrasolar planets, requires careful attention to 
various fine instrumental effects such as parasitic fringes. This constitutes an important motivation 
for optimizing the design of future planet-detecting interferometers like the future ground-based 
instrument of the VLTI, MATISSE, or the future NASA spatial nulling interferometer, FKSI. In 
fact the calibration of this parasitic effect by reference star seems very difficult at this level of 
required precision. Moreover the use of a phase abacus is not possible since the parasitic phase is 
an object-dependent quantity. Nevertheless in the marginal case of an unresolved source, it appears 
that the parasitized coherent flux can be calibrated by the parasitized photometry. Therefore we 
propose some solutions to prevent parasitic interference effects, and which should be taken into 
account during the design stage of a future instrument. 
We suggest the following design solutions : 

• To maintain the beams independence in certain parts of the instrument. For example, the 
left panel of Fig. 2 shows the crosstalk occuring if a spatial filtering plane is common to all 
beams. The solution would be to use an independent pinhole for each beam. 

• To separate the paths of each beam by a careful baffling everywhere it is possible inside the 
instrument. 

• To transport the beams in a non-co-phased way except just before the recombination. These 
fixed path differences maintain the beams out of the coherence length in order to prevent any 
parasitic interference when cross-talk is occuring. 

As a final conclusion, it appears that, up to now, little attention has been paid on the phenomenon 
of parasitic interference. This issue has been formalized here in a general multi-axial recombination 
scheme, with special attention on the differential phase. This work has b een conducted in the 



framework of the AMBER instrumen t (IPetrov and Amber Consortium! 120031) an d especially of its 



successor, the MATISSE instrument (jLopez. et al.1 (|2008l ); lLagarde. et al.1 (|2008l )). 



In the case of other instrumental configurations (co-axial scheme, nulling interferometry, ...), com- 
plementary work would need to be performed in order to define requirements adapted to a larger 
range of interferometers. 



We thank our colleague, Pierre Antonelli, for the attention brought to the parasitic light issue, 
and for the fruitfull discussion motivated by the MATISSE instrument design study. We also thank 
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Table 1: Values of various parameters discriminating each phase spectrum extracted from 



Barman et all (|200ll ). 





Type of star 


Tint 


Separation star/planet 


spectrum 1 


dM5 


500K 


0.05 AU 


spectrum 2 


dM5 


500K 


0.5 AU 


spectrum 3 


dM5 


1000K 


0.005 AU 


spectrum 4 


dM5 


1000K 


0.1 AU 


spectrum 5 


G2 


500K 


0.3 AU 


spectrum 6 


G2 


500K 


1 AU 


spectrum 7 


G2 


1000K 


0.05 AU 


spectrum 8 


G2 


1000K 


1 AU 



This preprint was prepared with the AAS IATgX macros v5.2. 
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Fig. 1. — Principle diagram representing the observation of an unresolved source by a 2-telescope interferometer. 
Fringes are obtained in the detector plane within a fixed fringe envelope corresponding to the point spread function 
of the telescopes. 8 represents the residual achromatic piston that will affect the fringes within the signal envelope. 
The area where crosstalk between beams occurs, thus causing a parasitic interference, is simply represented by a 
'module' located just before the beam recombiner. The fringe envelopes and patterns are normally overlapped at the 
same place on the detector but discriminated here for a matter of clarity. 
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Fig. 2. — Two examples of stray light and cross-talks between beams producing a parasitic interference : 
The left figure showing a beam overlapping due to diffusion effects, and the right one showing a cross-talk 
in a pupil plane due to the diffraction introduced by a common spatial filter. 
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Fig. 3. — Comparison between a$, the amplitude in L band of "J> po i nt (A) (noted 'intrinsic' in the figure), and a par *, 
the amplitude in L band of Xpoint(A) — $ P omt(A) (noted 'parasitic' in the figure), with respect to e" 2 . Each panel 
corresponds to a piston value (noted S in the paper) which is written in fraction of the central wavelength of the L 
band (A = 3.5/im). The A/500 and A/30 cases correspond approximately to the typical specifications that would be 
respectively achieved by a spatial and ground-based fringe tracking device (S w 5nm and lOOnm). 
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Fig. 4. — Comparison between the 'unpistonncd' intrinsic phase $ syst (A) (solid line) calculated from the number 
7 spectrum, and different 'unpistonned' parasitic phases Xsyst(A) — $ syst (A) calculated with different values of 8, 
expressed in fraction of A = 3.5/im. The A/500 and A/30 cases correspond approximately to the typical specifications 
that would be respectively achieved by a spatial and ground-based fringe tracking device (8* « 2nm and lOOnm). A 
parasitic flux factor of 10% has been considered here. 
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Fig. 5. — Amplitude, in L band, of the intrinsic phase of a stellar system with a hot Jupiter, a syst (horizontal 
solid line), and amplitude in L band of the related parasitic phase, a par (other dashed and dotted lines), plotted with 
respect to e" 2 . Different values of <5» ranging from A/500 to A/5 are considered, and each panel is related to a different 
hot Jupiter spectrum. The A/500 and A/30 cases correspond approximately to the typical specifications that would 
be respectively achieved by a spatial and a ground-based fringe tracking device (<5* « 2nm and lOOnm) . 



